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ABSTRACT 

We investigate forced liquid flows through open 
capillary channels with free surfaces experimentally. 
The experiments were performed under low grav- 
ity conditions in the Bremen Drop Tower and on 
board the sounding rocket TEXUS-37. Open cap- 
illary channels (vanes) are used in surface tension 
tanks to transport the propellant and to provide a 
flow path for the bubble-free liquid supply to the 
thrusters. Since the free surfaces can only withstand 
a certain pressure differential between the liquid and 
ambient, the flow rate in the channel is limited. The 
maximum flow rate is achieved when the surfaces 
collapse and gas is ingested into the outlet. Since 
experimental and theoretical data of this flow rate 
limitation is lacking, the safety factors for the ap- 
plication of vanes in surface tension tanks must be 
unnecessary high. 

The aim of the investigation is to determine 
the maximum liquid flow rate and the correspond- 
ing critical flow velocity. The characteristic non- 
dimensional parameters, Ohnesorge number, and 
gap ratio, cover a wide range of usual vanes. For the 
theoretical approach a one-dimensional momentum 
balance was set up. The numerical solution yields 
the maximum volume flux and the position of the 
free surface in good agreement with the experiments. 

INTRODUCTION 

Open capillary channels are used in space to 
transport and position liquids. An open capillary 
channel is a structure of arbitrary geometry which 
establishes a free surface flow when the hydrostatic 
pressure vanishes and surface tension forces become 
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dominant, e. g. within the compensated gravity en- 
vironment of an orbiting spacecraft. The free surface 
of the liquid flow is always exposed to the ambient 
pressure and the liquid is fixed inside the channel due 
to capillary forces. Prominent applications are heat 
pipes and surface tension tanks of satellites. In the 
latter case, open capillary channels, so called vanes, , 
provide a flow path from the bulk of propellant to 
the tank outlet or a refillable reservoir. Simple vanes 
consist of a thin sheet metal mounted parallel or per- 
pendicularly to the tank shell. Whenever the vane 
is in contact with the propellant the vane fills itself 
due to the capillary forces. Driven by a gradient of 
the capillary pressure between the liquid head in the 
vane and the bulk liquid, the propellant is conducted 
to the tank outlet. Thus, vanes are a simple device 
for the direct supply to the tank outlet wherever the 
bulk propellant is located in the tank. 

The major problem of the direct use of vanes 
is the limitation of the flow rate, which causes the 
liquid surface to collapse. Since the physics behind 
the effect of flow rate limitation in open capillary 
channels is not well-known, the design of vanes is 
based on estimations. The only work in this field 
are model computations by Jaekle 1 and Der 2 . To 
address the lack of experimental data we developed 
an experiment to investigate the limits of flow rates 
in open capillary channels. 

To identify the problem we consider an open 
channel as shown in Fig. 1. The channel consists of 
two parallel plates of distance a, width b and length 
l. The inlet and the outlet of the channel are each 
connected to a duct with a closed cross section. A 
constant volumetric flow Q in the direction of the 
longitudinal channel axis x is assumed. The liquid 
flow path is bounded by free surfaces between the 
edges of the plates. Both liquid surfaces are exposed 
to the ambient gas pressure p ai which is kept con- 
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Figure 1: Schematic drawing of the capillary channel 
consisting of two parallel plates. The volume flux Q 
through the channel along the flow path x leads to 
an increasing surface curvature. 


stant as well as the temperature. Depending on the 
applied volume flux, the liquid pressure decreases in 
the flow direction x due to flow losses. To achieve 
stationary flow conditions the difference between the 
liquid pressure p and the ambient pressure p a has 
to be balanced by the capillary pressure of the free 
surface. For that reason the surface curvature in- 
creases in the flow direction and the flow path be- 
comes smaller. Since the free surface can only with- 
stand a certain differential pressure, the flow rate in 
the channel is limited to a certain value. The max- 
imum rate is achieved when the surfaces collapse at 
the end of the capillary channel and gas ingestion 
occurs at the channel outlet. The aim of the investi- 
gations is to determine the maximum flow rate and 
the corresponding critical velocities. 

THEORETICAL APPROACH 

For the analysis we assume a one- dimensional 
flow along the channel axis x characterized by the 
mean velocity v and the liquid pressure p (see Fig. 
1). The origin of the coordinate system is located in 
the center of the cross section at the channel inlet. 
The basic equations to be solved are the momentum 
equation 

dp + pvdv -r du]{ — 0 (1) 

and the equation of continuity 

d{Av) = 0, (2) 


both given in differential form. The term dw( of 
Eq. (1) considers the pressure losses due to friction 
and the change of velocity profile in the entrance 
region, p denotes the density of the liquid and A is 
the cross section of the flow path perpendicular to 
the longitudinal axis x. We neglect the hydrostatic 
pressure since the experiments were operated in an 
environment of low gravity of 10~ 4 <? TO 5 g , with g 
representing the gravity acceleration on earth. 

At first a relation between the liquid pressure 
and the position of the free surface 2 is derived. The 
liquid pressure p is related to the mean curvature H 
of the free surface by the Gauss- Laplace equation 

p-p* = - 2(TH = - (T {-k + ~k)' (3) 

Herein o is the surface tension, R\ and R 2 are the 
principal radii of curvature. Since the ambient pres- 
sure p a is constant the pressure gradient becomes 

dp = —2(rdH. (4) 

In general, the mean curvature of the surface is cal- 
culated by 

x ( Z x + l) Z VV “ 22 x 2 y 2 xy + (I + Z 2 )z xx 

w(x ' y) (TTTT^ 

(5) 

(see e. g. Bronstein 3 ). 2 is the position of the sur- 
face defined by a function z = /(x,y), with x e [0,/] 
and y € [-a/2, a/2]. The indices in Eq. (5) denote 
partial derivatives: 2 * = dz/di , Zij = d 2 z/didj , with 
ij € {.r,y}. Eq. (5) can be simplified by applying 
the basic assumption of a one-dimensional flow that 
the liquid pressure variation is negligible over each 
cross section. With this assumption the surface cur- 
vature in the cross section is constant and defined 
by the radius of curvature R shown in Fig. 2. In 
that case the radius R and the cross section area A 
are calculated from the distance k(x) by geometrical 
relations 


R = 91W ( 6 ) 

A = Q2(R). (7) 


k is the innermost line of the surface 


k{x) = z(x,y)| y _o ( 8 ) 

in the plane of symmetry at y = 0, which was de- 
tected in the experiment. The simplification of Eq. 
(5) occurs since the derivative of 2 with respect to 
y vanishes along the innermost line 


dz 

dy 



(9) 
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Figure 2: Cross section of the capillary channel with 
radius of curvature /?, contour length fc, and sur- 
face deformation h. The observation is performed 
perpendicular to the plates. 


With Eq. (8) and Eq. (9) Eq. (5) reduces to 


pressure loss in a channel of two parallel plates pro- 
posed by Sparrow et. alA A: pe strongly depends 
of the position along the x-axis and vanishes if the 
entrance length is reached. 

Substituting Eq. (1) by Eq. (4), Eq. (12) and 
Eq. (13) the final differential equation for the change 
of pressure along the channel axis 

o dH &L d A _ pQ 2 ( fc Pf + fc P e) = „ lu) 

2a dx A 3 dx 2 D h Re Dl , A 2 ( 

is obtained. To get the dimensionless form the 
lengths are scaled by the hydraulic diameter D h and 
the channel width b. The liquid vel ocity is scaled by 
the characteristic velocity t’o = 2y/crJpD\ h, which is 
known from the capillary imbibition of liquid in a 
channel of two parallel plates. As shown by Dreyer 6 
this velocity cannot be exceeded. With x = 4x/Dh, 
r = 4//D h , 4 H = D h H, k = 4 k/D h , R = 4 R/ D h , 
X - 2 A/D h b, V = v/vo, Q = Q/ \JvD h b 2 /p the 
dimensionless form of Eq. (14) reads 


2H(x y y = 0) = 


+ 


(1 + kl)i i?(l + fc 2 )i 


( 10 ) 


With the definition of the volume flux Q from Eq. 

( 2 ) 

Q = Av (11) 

the second term of Eq. (1) can be rewritten as 


, Q 2 , A 

pv dv = —p—aA. 


( 12 ) 


The pressure loss dw[ in Eq. (1) consists of 
two parts, the laminar viscous pressure loss and the 
additional entrance pressure loss, which is a result 
of the change of the velocity profile from the en- 
trance profile to the parabolic velocity distribution 
(PoiSEUlLLE flow). The pressure loss we apply is 


dwt = 


pQ 2 


2 D h Rem A 2 


[k pf + k pe ]dx. (13) 


Neglecting the factor k pe , Eq. (13) yields the laminar 
pressure loss for two parallel plates obtained from 
an exact solution of the Navier-Stokes equation 
(White 4 ). Therein k pf = 96 is the friction factor, 
Du = 2 a the hydraulic diameter, and Reo h = Duv/v 
is the Reynolds number. The kinematic viscosity 
is denoted with v. The additional pressure loss due 
to the transformation of the entrance velocity pro- 
file to the laminar parabolic velocity distribution is 
taken into account by the factor k pe . This factor, 
which is not discussed in detail in this paper, was 
derived from an analytical relation of the entrance 


dH' Q ' 2 dX 

dx' + A * 3 dx * 


-^Q'(V+g=0. (15) 


The dimensionless parameters are the Ohnesorge 
number Oh = p/y/Dupv and the aspect ratio A = 
b/a occurring in the non-dimensional relation of the 
cross section A = g^{R )• Eq. (15) has to be solved 
numerically taking into account the boundary con- 
ditions which are 


= *1 = A, (16) 

\r= o 

= 2 Hi (17) 

the position of the surface at the inlet and outlet of 
the channel Eq. (16) and the surface curvature at 
the channel inlet Eq. (17). The numerical solution 
yields the position of the liquid surface 

k =/(A,Oh, l\Q m ). (18) 



Eq. (18) shows that the surface position depends on 
the adjusted volume flux Q . For each set of A, Oh, 
and r several model computations with different Q 
are performed. We assume that the maximal vol- 
ume flux is reached when the numerical integration 
diverges. As the maximal theoretical volume flux 
Q'rit we define the maximal volume flux Q lead- 
ing to the convergence of the numerical algorithm. 
Obviously the following relations 

Q'cnL = /i(A,Oh,0 (19) 

v; rit = /a(A,Oh,0 (20) 
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are valid. i** rit is the corresponding critical mean ve- 
locity u‘ rit = Q’ rit /i4^ in which occurs at the smallest 
cross section A mm of the liquid flow. 

From another point of view we expect that the 
limitation of the flow rate occurs due to choking. 
The theory of choked flow predicts a limiting veloc- 
ity corresponding to a characteristic signal velocity 
of the flow. In compressible gas duct flow's the lim- 
iting speed is defined by the speed of sound. The 
characteristic number is the Mach number (Ma), 
and the maximal flow passes through the duct when 
Ma = 1 is reached locally. 7 In open channel flow 
dominated by normal gravity the speed of shallow 
water waves correspond to the limiting velocity, and 
choking occurs when the Froude number (Fr) be- 
comes unity. Both wave speeds are derived from the 
general form 

c= s/Adpjd{pA). (21) 


Eq. (21) is the wave speed of longitudinal waves 
in compressible flow with variable cross section. 8 A 
close analogy exists between both flow types which 
might be expandable to open capillary flows investi- 
gated in the present work. To show that the flow rate 
limitation is caused by choking, we scale the mean 
liquid velocity by the characteristic wave speed of 
the open capillary channel. Applying Eq. (3) and 
Eq. (10) to Eq. (21) the wave speed of the capillary 
open channel is 


v c 


a A dR 
~p~W~dA' 


( 22 ) 


Here we neglected the surface curvature in flow di- 
rection assuming 


k xx - 0, k x « 1, (23) 


which is valid in pieces along the channel in sufficient 
distance to the inlet and outlet. With the relations 
Eq. (6) and Eq. (7) the wave speed can be calculated 
from the experimental determined surface contour 
k(x). Scaling the mean liquid velocity v = Q/A by 
Eq. (22) the dimensionless parameter 

We = Q/Av c = v/v c " (24) 

is defined, which is a Weber number. Due to the 
above mentioned analogy we expect that the We- 
ber number tends towards unity when the flow rate 
becomes maximal. 


EXPERIMENTS 

The objective of the experimental investigations 
Was to determine the maximum liquid volume flux 


Degassing valve 



Figure 3: Schematic drawing of the experiment 

setup for the drop tower. During the microgravity 
period the capillary channel fills itself and pumping 
starts when the meniscus reaches the screen. 
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which can be sucked through open capillary chan- 
nels, consisting of two parallel plates. To eliminate 
the gravity effects and to establish a liquid path be- 
tween the parallel plates by means of the surface 
tension only, the experiments were performed un- 
der microgravity conditions in a drop tower and in 
a sounding rocket. In addition several preliminary 
investigations and component tests for the sound- 
ing rocket experiment were carried out in two flight 
campaigns in a parabolic aircraft. 

The first experiments were performed in a drop 
tower and based on the experiments of capillary liq- 
uid rise between two parallel plates, Dreyer et. al. 6 9 . 

In essence the experiment setup consisted of a circu- 
lar container, which was filled with a test liquid up 
to a certain height as shown in Fig. (3). Inside the 
container a channel of two parallel plates was fixed 
on the bottom, so that under normal gravity condi- 
tions the longitudinal channel axis was perpendicu- 
lar to the liquid surface. On the top of the channel 
a withdrawal device consisting of a piston with step- 
per motor was mounted. The experiment setup was 
integrated in a drop capsule at the Bremen Drop 
Tower providing an experiment time of 4.74 s at a 
residual acceleration of < 10 _5 g. Before the drop the 
withdrawal device was filled with the test liquid up 
to the screen, but the gap between the plates upside 
the liquid surface was empty. With the release of the 
capsule the liquid started to rise between the plates. 
This liquid movement is caused by the capillary pres- 
sure gradient between the liquid surface inside the 
container and the liquid surface inside the gap of 
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the plates. When the liquid meniscus contacted the 
withdrawal device the stepper motor started and a 
defined constant volume flux was withdrawn from 
the channel by means of the piston. Depending on 
the adjusted volume flux two typical observations 
were made. Below the critical value of the volume 
flux Qcrit (Q < Qcrit) the flow between the plates is 
stationary. Above the critical value (Q > Qcrit) the 
surface collapses and gas ingestion occurs. Based on 
this observation* we define the maximum volume 
flux as 

Qcrit = 0.5(Qbc + Qac) » (25) 

where Qb c is the maximum experimental volume flux 
leading to a stable flow and Q ac the minimum exper- 
imental volume flux leading to a the surface collapse. 

Several experiments with various liquids and dif- 
ferent channel properties were performed. Since the 
result of this work shall be useful for surface ten- 
sion tanks, the experimental parameter were chosen 
with respect to propellant data (N2O4 and MMH) 
and tank geometries. Satellite tanks usually oper- 
ate in the temperature range T = 0 . . . 40 °C and 
typical geometries are: gap distances a « 1 ... 10 
mm, plate breadths 6 « 10 ... 30 mm. This leads 
to the parameter field given by the combination 
0.427 • 10~ 3 < Oh < 5.881 • 10” 3 and 3 < A < 10. 
The experiment parameters are given in Table 1 and 


Table 1: Experiment parameters for Bremen drop 
tower experiments with experiment time of 4.74 s. 


# 

Fluid 

a/cm 

6/cm 

//cm 

1 

SF 1.0 

0.5 

3 

9.5 ± 0.05 

29 

FC-77 

0.2 

2 

7.5 ±0.1 

32 

SF 0.65 

0.3 

1 

9.45 ± 0.05 

33 

SF 0.65 

0.3 

1.5 

9.6 ±0.05 

34 

SF 0.65 

0.3 

3 

9.45 ± 0.05 

35 

SF 0.65 1 

0.5 

2.5 

9.5 ± 0.05 

35a 

SF 0.65 

0.5 

2.5 

4.8 ±0.05 

36 

FC-77 

0.2 

1 

7.5 ±0.1 


Table 2. All fluids show perfectly wetting behavior 
at the capillary walls. Laminar flow conditions were 
expected for all experiments since the REYNOLDS 
number was sufficiently small. The combination of 
flow length l and maximum volume flux Q must be 
chosen with respect to the experiment constraints: 
overall experiment time (4.74 s), capillary rise time 
to fill the channel compared to pumping time, op- 
tical resolution of the surface contour (depends on 

•A detailed explanation is given below, when the sound- 
ing rocket experiment is discussed. 


Table 2: Dimensionless parameters for the drop 

tower experiments. 


# 

A 

Oh 

— ■ 

1 


1 

5.99 

0.00238 

0.045 

534 

29 

3.33 

0.0047 

0.177 

156 

32 

3.33 

0.00197 

0.062 

578 

33 

5 

0.00197 

0.063 

600 

34 

10 

0.00197 

0.062 

625 

35 

5 

0.00152 

0.029 

909 

35a 

5 

0.00152 

0.015 

977 

36 

5 

0.00472 

0.177 

151 


the gap distance a). Due to the short experimental 
time the Ohnesorge number was restricted to Oh 
> 1.52- 10“ 3 (#35). 

The maximum volume flux was determined for 
each parameter set. Theoretical model computa- 
tions were performed under the assumption of Eq. 
(23) and a different factor k^. The calculations show 
a good agreement with the experiments with respect 
to the maximum volume flux and the predicted sur- 
face contour. However, due to the transition from 
the self-driven capillary flow to the forced flow iner- 
tia effects occurred which made an exact approach 
of the critical value impossible. Furthermore, the es- 
tablished volume flux had to be kept constant since 
the short experimental time did not allow an increase 
of the volume flux. More details regarding the exper- 
iment setup and the results are discussed by Dreyer 
et. al. 10 . 

To operate under longer /ig-time an experiment 
module was flown on the sounding rocket TEXUS- 
37 providing a microgravity environment of 10~ 4 g 
in all axis for approximately 6 minutes. The module 
was equipped with a pump system that allowed the 
increase of the volume flux in small steps (quasistatic 
approach) to minimize the transient inertia effects. 

A schematic drawing of the experiment cell is 
shown in Fig. 4. It consists of the liquid reservoir 
with compensation tube and the observation cham- 
ber with the channel (a = 5 mm, b = 25 mm, l — 47 
mm). The flow is established by two gear pumps. 
One pump supplies the reservoir through a circular 
gap on the bottom of the reservoir (volume flux Q 1). 
Via the nozzle, the liquid is conducted into the cap- 
illary channel. At the channel outlet the volume flux 
Q 2 is withdrawn by the suction pump. The differ- 
ence of both volume fluxes caused by fluctuations 
of the rotation speed and varying liquid slip inside 
the pump is compensated by the compensation tube. 
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Figure 4: Schematic drawing of the TEXUS experi- 
ment module. 

Furthermore the liquid meniscus in the tube sets the 
pressure in the reservoir and defines the boundary 
condition at the channel inlet Eq. (17) needed to 
solve Eq. (15). To prevent the meniscus disappear- 
ing into the reservoir the supply pump was adjusted 
to pump a 3% higher volume flux than the suction 
pump. With this adjustment, the liquid rose con- 
tinuously in the compensation tube. Each pump is 
connected to a separate reservoir for the test liquid 
supply and the storage of withdrawn liquid. 

The experiment was observed by three CCD 
cameras. The video signals were down linked to the 
ground and recorded on a VCR (S-VHS quality). 
Additionally the camera view of the compensation 
tube was recorded on board by a DV camcorder. 
Two cameras were used for the observation of the 
channel (made from quartz glass) yielding a resolu- 
tion of approximately 0.05 mm/pixel. Both sections 
of image had an overlap of 10 mm. On the front plate 
markings were etched for the calibration of pictures 
and the check of the optical resolution. 

The experiment was controlled manually dur- 
ing the flight of the sounding rocket. All hand- 
controlled procedures necessary for the TEXUS ex- 
periment control had been practiced during two par- 
abolic flight campaigns. In the initial configuration 
before the lift off the compensation tube and the noz- 
zle were each closed with a gate valve, both pumps 
were off. Once under microgravity conditions, the 


capillary channel was filled. To avoid wetting out- 
side the channel caused by an uncontrollable liq- 
uid jet, the channel was not filled directly by the 
pressure pump but by the self-driven capillary rise, 
which was well-known from the drop tower experi- 
ments. For that reason, the compensation tube was 
filled first. After opening the gate valve at the noz- 
zle the liquid rise occurred and the channel filled ■ 
itself. When the liquid reached the channel outlet 
it was withdrawn by the second pump. By starting 
the first pump a continuous flow through the open 
channel was established. 

The filling procedure took approximately 30 s. 

In the remaining experimental time the volume flux 
was increased up to the critical value. In contrast 
to the drop tower experiment the change of volume 
flux by each step is small, and the inertia effects are 
sufficiently smaller since the channel is supplied by 
the two synchronized pumps. The experiment pro- 
vided the options to change the volume flux in 0.1 
ml/s or 0.05 ml/s steps. At first the bigger step size 
for a rough approach was chosen to save experiment ^ 
time. After the gas ingestion was observed the vol- 
ume flux was decreased and a second approach was 
performed at small steps of volume flux. The typical 
observations, which are the same for the drop tower 
experiments, are shown in Fig. 5 and Fig. 6. The 
camera view is perpendicular to the top plate. Since 
the plates are transparent the free surface appears as 
a dark contour, which corresponds to the distance h 
defined in Fig. 2. Below the critical value of volume 
flux Qcrit {Q < Qcrit) the flow between the plates is 
stationary. The free surface is stable and the corre- 
sponding innermost line of the surface is constant in 
time, k = k{x). For Q > Q c rit the surface collapses 
and gas ingestion occurs. The free surface becomes 
time-dependent, k = k(x, £)• The withdrawn total 
flux now is a superposition of the maximum liquid 
volume flux and an additional volume flux of gas. 

Prior to the experiment, we performed three- 
dimensional model computations using the finite el- 
ement code FIDAP 11 . The computations were nec- 
essary to optimize the flow conditions inside the liq- 
uid reservoir and the nozzle as well as to calculate 
the pressure loss caused by both components. With 
the knowledge of the pressure loss the pressure at 
the channel inlet and the related boundary condition 
from Eq. (17) (necessary to solve Eq. (15)) is de- 
fined. Both liquid surfaces in the compensation tube 
and in the channel are exposed to ambient pressure. 
For this reason the inlet pressure is determined by 
the ambient pressure minus the static capillary pres- 
sure of the liquid surface in the compensation tube 
and the dynamic flow losses, which 


6 




nnnnnnpn 

i.ik ii i i 



Q = 6 ml/s <3 = 7 ml/s Q = 8 ml/s Q = 9 ml/s 


Figure 5: Stable volume flux during the TEXUS experiment, Q < Q C rit> The camera view is perpendicular 
to the front plate as shown in Fig. 2. Since the plates are transparent the free surface appears as a dark 
contour. The flow direction is from the bottom to the top. With increasing volume flux the surface curvature 
increases along the channel axis and the flow path becomes smaller. 



t = 197.02 s t= 197.46 s t = 197.86 s t = 197.94 s 


Figure 6: Unstable liquid flow in a capillary channel between two parallel plates during the TEXUS exper- 
iment. The adjusted volume flux lies above the critical value Q > Qcrit • The differential pressure between 
the liquid and ambient cannot be compensated by the capillary pressure of the liquid surface. The surface 
collapses and gas ingestion occurs. 
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was fitted to 


60 


2H 0 = k\ + k 2 Q + ^’3 Q “ (26) 

in terms of the surface curvature. The constant 
k x = 0.222 considers the capillary pressure of the 
compensation tube (inner diameter 45 mm). Since 
the test liquid is perfectly wetting, the shape of the 
surface inside the tube is a section of a sphere cu- 
lotte. The factors k 2 = 0.128 and k 3 = 0.725 repre- 
sents the dynamic losses due to friction and convec- 
tive acceleration. 

RESULTS 

During the 6 minute flight the experiment mod- 
ule worked well. The filling procedure and all ex- 
periment procedures were performed on schedule. 
No leakage of liquid into the observation cham- 
ber occurred, even though the pumps were stopped 
abruptly. After the flow through the channel was 
established, the volume flux was increased up to the 
critical value twice. With the first approach at the 
bigger step size the surface collapsed at Q = 9.1 
ml/s. At the same value of Q the second approach 
with smaller steps yielded a stable flow, which could 
not be increased further since the experimental time 
expired, thus the critical volume flux by the defin- 
ition of Eq. (25) was not obtained. However, this 
effect shows how sensitive the flow responds to iner- 
tia effects at this stage. 

Several stable flow configurations at different vol- 
ume flux Q below the critical value were observed 
as shown in Fig. 5. As expected the distance h 
increased with increasing volume flux and the liq- 
uid flow path became smaller. The highest volume 
flux leading to a stationary flow could be realized 
at Qbc = 9.1 ml/s. The recorded video pictures are 
processed by digital imaging tools so that the inner- 
most line k(x) is available for each stable volume flux 
(Q < Qcrit)« For each Q during the experiment Eq. 
(15) with the boundary conditions from Eq. (26) 
was solved. The comparison with the surface con- 
tour data is shown as an example in Fig. 3 for a low 
and a high volume flux at Q — 6 ml/s and Q = 9 
ml/s, respectively. The evaluated surface positions 
at different times are represented by the black ter- 
raced lines. Slight surface vibration caused by the 
gear pumps are observed. The solid line is the the- 
oretical prediction of k. The figures show that com- 
putations predict the volume flux and corresponding 
the surface position in good agreement with the ex- 
perimental data. 

To show the influence of the Weber number 
from Eq. (24), the waves speed v c Eq. (22) and the 




Table 3: Evaluated surface contour from the TEXUS 
experiment (thin lines). Slight surface vibrations 
caused by the gear pumps are observed. The thick 
solid line shows the theoretical prediction (a: Q = 6 
ml/s, b: Q = 9 ml/s). 
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Figure 7: WEBER number as function of the ad- 
justed volume flux. At We = 1 the flow is choked 
and a further increase of volume flux is impossible. 

mean liquid velocity v = Q/A are calculated from 
the experimental data of the surface position k(z). 
Using the relations Eq. (6) and Eq. (7) both veloci- 
ties are determined. Since the effect of choking is ex- 
pected to appear at the location of minimum liquid 
pressure, the WEBER number was calculated at the 
smallest cross section A m , n = A(/c m j n ). The analy- 
sis shows that with increasing volume flux the liquid 
velocity increases, but the characteristic wave speed 
decreases. Therefore, the Weber number tends to- 
wards unity as the theory of choking predicts (Fig. 
7). At the maximum stable volume flux (Q^ c = 9.1 
ml/s) the corresponding Weber number yields We 
= 0.86. Fig. 7 indicates that the flow runs very 
close to the critical condition, since a small change 
of volume flux causes a large variation of the Weber 
number at this stage. 

OUTLOOK 

As we have demonstrated in the drop tower and 
sounding rocket experiments, much longer micro- 
gravity durations are necessary to fully span the pa- 
rameters discussed in this paper. Furthermore, very 
small incremental increases in flow rates are required 
to minimize the transient inertial effects. In order to 
meet these requirements, this experiment has been 
tentatively selected as a definition-phase flight ex- 
periment for the ISS. The experimental hardware 
will be developed and built at the Center of Applied 
Space Technology and Microgravity (ZARM) and in- 
tegrated into the Fluids Integrated Rack (FIR). The 
FIR is a component of the Fluids and Combustion 


Facility (FCF) and will be located in the US labo- 
ratory module, Destiny. 

We will take advantage of both the 30 fps high 
resolution digital black and white camera package 
and the Ultra High Frame Rate (UHFR) camera. 
We will also make use of the Command and Data 
Management Subsystem (CDMS) to control and 
measure liquid and gas flow rates, temperatures and 
pressures. Integration and testing into the FIR is 
tentatively scheduled for 2006. 

CONCLUSION 

The aim of the experiments was to investigate the 
flow rate limitation in open capillary channels. In 
spite of the applications of capillary channel (vanes) 
in surface tension tank technology the basic physics 
leading to the flow rate limitation is not well-known. 
Experimental and theoretical data on this research 
field are rare. To eliminate the hydrostatic pressure 
the experiments were performed under the micro- 
gravity environment of the Bremen Drop Tower and 
on board the sounding rocket TEXUS-37. The inves- * 
tigated dimensionless characteristic numbers cover a 
wide range of usual vanes. For the theoretical ap- 
proach a one-dimensional momentum balance was 
set up, which yields the maximal volume flux and the 
position of the free surface k(x). Both the calculated 
volume flux and the surface position are in good 
agreement with the experimental data. Furthermore 
the experiment on board TEXUS-37 showed that 
the limitation of flow rate occurs due to choking 
indicated by the Weber number. Since the char- 
acteristic wave speed v c was derived from the gen- 
eral waves speed Eq. (21) the capillary open chan- 
nel flow is analogous to Mach number and Froude 
number problems. For the entire understanding of 
the mechanism of flow rate limitation in open capil- 
lary channels the parameter field to be investigated 
must be expanded. To keep the influence of inertia 
effects smaller the increments of volume flux have to 
be reduced leading to longer experimental times. In 
order to meet these requirements, this experiment 
has been tentatively selected as a definition- phase 
flight experiment for the ISS. 
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NOMENCLATURE 
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a distance between the parallel plates 

b breadth of parallel plates 

g gravitational acceleration 

c longitudinal wave speed in fluids with vari- 

able cross section 
h height of circle segment 

k position of the innermost surface line 

A: min minimum position of the innermost surface 
line 

k v { factor of the laminar pressure loss, k p f = 96 

k pe factor of the entrance pressure loss 

l length of capillary channel 

p pressure 

p a ambient pressure 

v average velocity in direction of the flow path 

vo characteristic velocity, t*o = \jlajpa 

v c longitudinal wave speed in open capillary 

flow 

v crit critical flow velocity 

W[ energy loss of the flow 

x channel co-ordinate, flow path co ordinate 

y channel co-ordinate 

z surface position 

A cross section area of the flow path 

A m j n minimum cross section area of the flow path 

Dh hydraulic diameter, Dh = 2a 

Fr FROUDE number 
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Oh OhnesorGE number, Oh= pjsJDhPv 

R radius of curvature 

Ri <2 main radii of curvature 

Reo h Reynolds number based on the hydraulic di- 
ameter, Reo h = Dhv/v 
Ma MACH number 

We Weber number, We— vj v c 

p density 

p dynamic viscosity 

v kinematic viscosity 

cr surface tension 

A gap ratio, A = bja 

* dimensionless quantiy 
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